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Abstract. - We use the microscopic weak coupling theory to predict the pairing state in 
superconductors of cubic, hexagonal, or tetragonal symmetry, where the order parameter is 
multicomponent, i.e., transforms according to either a 2-dimensional or a 3-dimensional repre- 
sentation of the crystal point group. We show that the superconducting phase usually breaks 
the time-reversal symmetry for singlet multicomponent superconductors. The superconducting 
order parameter for triplet superconductors in most cases turns out to be non-magnetic. 



While most superconductors are described by the Bardcen-Cooper-Schrieffer (BCS) weak 
coupling theory, which gives s-wave symmetry of superconducting state, in many other materi- 
als, such as ^He [1], UBeia, UPts [2], the high-Tc cuprates, Sr2Ru04 [3-5], and PrOs4Sbi2 [6], 
the pairing state is known to be unconventional. (See Refs [7,8] for a review). If the pairing 
state in an unconventional superconductor or superfluid transforms according to a multi- 
dimensional representation of the group, the order parameter is also multicomponent. For 
example, the p-wave order parameter in "^He is a 3 x 3 complex matrix, i.e., has 9 complex 
components. As a result, many different phases can arise in the superfluid state, such as the 
A-phase [1], which is not invariant with respect to rotations or time reversal operation, and 
the B phase, which is both rotationally and time-reversal invariant. Unconventional super- 
conductors usually have a gapless excitation spectrum, which has important consequences for 
thermodynamics. Instead of describing the fascinating properties of unconventional supercon- 
ductors, we refer the reader to some excellent books and reviews on the subject [2,5,7-10]. 

Although the details of the mechanism are essential for strong coupling Eliashberg-type 
calculations, the knowledge of the energy spectrum and the coupling constants for quasipar- 
ticles mediating superconductivity is not required [11] in the weak coupling approach, which 
is adopted below. A microscopic description of unconventional superconductors can be useful 
for analysis of experiments in new superconducting materials, where unconventional pairing 
mechanism is suspected. The usual suspects are materials with strong Coulomb correlations, 
such as the heavy fermions or high-Tc cuprates. In such materials conventional phonon mech- 
anism gets suppressed [12]. 

The microscopic analysis is somewhat less general then the Ginzburg-Landau (GL) ap- 
proach of Refs [7, 8] . For example, weak coupling BCS-type theory always gives the B-phase 
of superfluid '^He [13]. The non-BCS spin fluctuation feedback effect, when the interaction 
itself depends upon the superconducting ground state, is required to stabilize the A phase [14]. 
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Most, if not all, unconventional superconductors are characterized by the presence of strong 
Coulomb correlations, so the applicability of weak coupling approach for such materials is 
questionable. Nevertheless, this approach is often used theoretically as a starting point. 

In what follows we analyze multicomponent superconductors, i.e., superconductors belong- 
ing to a 2-dimensional or 3-dimensional representations of the cubic, hexagonal, and tetragonal 
crystal point groups, where the superconducting phase is not uniquely determined by the ir- 
reducible representation [7,8]. Examples of superconductors where multicomponent state is 
suspected are, among others, UPta [2] (hexagonal symmetry), Sr2Ru04 [3, 4] (tetragonal sym- 
metry), and, more recently, PrOs4Sbi2 [6](cubic symmetry) and Na2:Co02 • 2/H2O [15,16]. 
Like in superfluid '^He, many phases can possibly appear in these superconductors, according 
to the GL theory [7,8]. However, only one of these phases is realized in the microscopic weak 
coupling theory. The analysis done below is well known [10,11], and has been applied to many 
particular cases, such as the d-wave superconductivity in high-T^-s [17], UPts [18], or general 
multi-band phonon mechanism [19]. 

As usual [20], we assume that electrons have a weak short-range attractive interaction, 
U{r,r'). Since the phonon cutoff frequency is much less than the Fermi energy, ujd ^ ep, 
only the quasiparticles in the vicinity of the Fermi surface participate. Then the interaction 
can be expanded in a complete set of basis functions, xim{p), where index / corresponds to 
different irreducible representations of the point group, while index m enumerates the basis 
functions of a given irreducible representation I. This expansion is similar to expansion in 
spherical functions for an isotropic model. The interaction in momentum space then takes the 
following form: 

u{p,p') = J2 UixiMximip'). (1) 

Im, 

where p, p' are on the Fermi surface. The gap equation for the superconducting order pa- 
rameter takes the following form: 

Ao^,™(p) = \Ui\xira{p) j -^Xlmip') ^ ^^„;3 (p' ; , (2) 

where Fap{p'] iw„) is the Fourier component of anomalous Gor'kov function, 

F^0{r-r',T-T') - -(T, (*„(r,r)*0(r',r'))), (3) 

and Ui < is a constant in Eq. |^ corresponding to a selected pairing channel. The two fermion 
operators inside brackets in Eq.Q anticommute. In the presence of the center of inversion, the 
behavior of xi{p) with respect to inversion, p — > p, determines the symmetry of Aq^(p) - 
even for singlet, and odd for triplet pairing [7,8]. The order parameter below Tc has the form: 

Aa/3,m(p) = l^{T){iay)aiiXlm{p), = 0, (4) 

Ac/3,m(p) = A(T)(cr„^ ■d)x(™(p), S^l. (5) 

Different phases for the multicomponent pairing state were determined from the GL analysis 
in Ref. [7,8]. The weak coupling approach completely determines the coefficients [11] in the 
GL functional. The integral gap equation Eq.j^l) written in the vicinity of Tc becomes the 
GL equation. In what follows we analyze the GL coefficients in degenerate representations 
of cubic, hexagonal, and tetragonal groups to determine possible phases that can arise in the 
microscopic theory. 
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The cubic group Oh - two-dimensional representations. - The order parameter for S=0 
and S=l can be written as 

where the basis functions are 

'f^^^k) oc kl + ekl + e^kl = ^ = 0, (7) 

ip^^^k) xk^ + eyky + e^zk^, 0'^^^ (k) = ip^^\k)* ; S = 1. (8) 

Here e = exp (27rz/3), and the wave functions are norniahzed to 1: 

(l^(^'(fc)P)F.5.-i, \m\' + \V2\' = l- (9) 
For singlet order parameter, we find: 



F = NnA^ In 



T 



Here A^o is the density of states for one spin direction, a is a coefficient determined from the 
Fermi surface averages: 

_{\0W{k)\^)p.s. ^^^^ 



(|(^(l)(fc)|2\2 



F.S. 

For the isotropic spherical shape of the Fermi surface we find that a = 10/7. While a will 
be different for an arbitrary Fermi surface, the form of the functional Eg. 1)10(1 is completely 
determined by the cubic symmetry, since certain averages, such as {ip^^'^ {k)'^) p s ^ must vanish 
by symmetry (in this case, rotation around the 3-fold axis). This functional leads to the phase 
(1, 0) with a magnetic d-wave order parameter, 

A{k) ockl + ekl + e^kl (12) 

which breaks T — > — T symmetry, and has point nodes along the diagonals of the cube. Note 
that in this case higher order terms are not necessary to determine the phase. In the triplet 
case, we obtain, for the spherical Fermi surface: 

j^(4)_(6) '''C(3)iVo^'' /, |4 , I |4 , I |2| |2^ ^'^C{5)NqA^ 3 *3 , *3 3\ /i q\ 

10^2^2 (1^1 1 +l^2| \V2\)- 2240(7rTe)4 ^"^'"^^ + V2) (13) 

The 6-order term is necessary in this case to distinguish between the two non-magnetic su- 
perconducting phases. The magnetic class turns out to be the trivial D4 ® R, with the order 
parameter 

d{k) (X 2xkx — yky — zk^. (14) 

This order parameter is nodeless, i.e., thermodynamic properties of this state are exponential 
at low temperature, similar to BCS s-wave. For a general Fermi surface we found that the 
sign of the 6-th order term is always negative. However, the coefficient k in the 4-th order 

oc |ryi|4 + 1,^2 14 + ^1^1 (15) 

depends on the shape of the Fermi surface. We find that —0.8 < k < 4, which makes magnetic 
class 0(02) {k > 2) also a possibility. 
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Three-dimensional representations. - For the four 3D vector representations, the basis 
functions can be chosen in the following form: 



FluiS = 1) 

F2g{S = 0) 
F2g{S - 1) 



kykz (ky k^)^ kzkx{k^ ^x)^ ^-cikx ^y)'i 



ykz zky-f zkx 



hkz 



ykx 



The Free energy can be written in the universal form: 



F = NqA^ In 



(rf-Tf^) + /3iNQA\Tf- rf^y + /?2 7VoA*|rf • rf|^ + P3N0A 



a) Fig. Wc find the following GL coefficients in Eg. 1)20(1 for a spherical Fermi surface: 

2205C(3) 



Pi 



35367r2T2 



132 



Pi 

2 



22 



Pi, 



(16) 
(17) 
(18) 
(19) 



(20) 



(21) 



which gives magnetic superconducting class D3{E), with a phase (l,e, e^), given in terms of 
the above basis functions. This phase is ferromagnetic, with orbital magnetic moment pointing 
along the 3-fold axis of the cube. The order parameter has 2 point nodes, at the points of 
intersection of the spontaneous 3-fold anisotropy axis with the Fermi surface. For a general 
Fermi surface, the phase is not completely determined. We found that 



^ -2' W,^-2' 



(22) 



This gives a line on more general phase space for the GL theory. Thus, magnetic phases 
D3{E), Di{E)^ and a non-magnetic phase Di{Ci) ® R are possible depending on the shape 
of the Fermi surface. 

b) F2g. For a spherical Fermi surface, we find: 



Pi 



5C(3) 
87r2T2 ' 



P2 



Pi 
2 



/33 = 0. 



(23) 



This case falls on the boundary between phases (l,e, e^) and (l,i,0), with accidental degen- 
eracy between D4^{E) and D3{E). For an arbitrary Fermi surface, we find: 



P2 
Pi 



P2>0, I > 
Pi 



5 

'4 



(24) 



(2) 

and the magnetic phases D3{E), D4{E), and non-magnetic {D2) ® R are possible. 

c) Flu and F2u- For both triplet representations we find, in case of spherical Fermi surface: 



Pi 



63C(3) 
807r2T2 



P2 



Pi 

3 



/33 = 0, 



(25) 



which falls on the boundary between 2 non-magnetic phases: (1, 1, 1) (class I?3(C3) ® R for 
Flu, D3(»R for i^2u) and (1,0, 0) (class D4{C4) ® R for class ^ {D2) (E) R for F2u)- For 
the general Fermi surface one of these non-magnetic phases will be realized. We find that 



/?3 ^^P2 
Pi ^ Pi 



0, /3i > 0, ^ < 0. 

Pi 



(26) 
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The hexagonal group -Dgft- ~ Multicomponent order is possible in hexagonal group, if the 
pairing wave function transforms according to one of the 2D irreducible representations. For 
El, the GL functional is given by Ea. (|20() . with the order parameter now a 2D vector: 

Eig{S = 0) : i>{k) oc rj^k^k^ + rjyk^ky (27) 
Eiu{S = 1) : d{k) cx ri^zk^ + rjyzky (28) 

For an arbitrary Fermi surface we find: 

^ = ^ , /3i > 0, /?3 = (29) 
Pi ^ 

This gives rise to the phase in both singlet and triplet cases, or the symmetry class 

Dq(E), which breaks the time- reversal symmetry, and possesses ferromagnetism. For the 
particular case of a spherical or cylindrical Fermi surface we find: 

ft(,..0^|||. AW.)^^^ (30) 

For E2g the GL functional is the same as for the 2D representation Eg in the cubic lattice. 
The basis functions are, of course, not the same: 

E2giS = 0): ^^^\k) {k, + ikyf , = (31) 

For an arbitrary Fermi surface the form of the functional, Ea. p()|l . stays the same, while the 
coefficient a depends on Fermi surface averages of the basis function. For example, a for the 
particular cases of spherical or cylindrical Fermi surface given by: 

a{sph) = — , a{cyl) = 1. (32) 

The resulting phase ip(k) cx (fc^ + iky)"^ belongs to the symmetry class Dg{C2), and is an 
orbital ferromagnet. Finally, for E2u triplet representation, with basis wave functions 

E2u{S^l): ip^'Hk)o,ix + iy){k,+tky), (^(^^ (fc) - (fc)* (33) 

the free energy functional is given by: 



F = In 



T 



NoA^Qr^if + \7j2f ) + aNoA*^^ (Iryi^ + 1^21"), (34) 



which gives either the p-wave phase xkx—yky (symmetry class D2®R), or xky+yk^ (symmetry 
class 02(02) (8> R)- Unfortunately, BCS does not distinguish between the two phases in the 
6-th order. Here a > is a coefficient, which depends on the shape of the Fermi surface. For 
a cylindrical and spherical Fermi surface, 

a{sph) — ^ , a{cyl) — 1. (35) 

The tetragonal group ■ ■ ~ Since the basis functions for the 2D representations E^ and 
Eg of the tetragonal group are exactly the same as for hexagonal and Eig representations, 
not surprisingly, the GL functionals and the resulting phases also turn out to be identical. 
The class in tetragonal group is D^lE), which allows ferromagnetism along z-axis. There are 
2 point nodes for the singlet, and a line of nodes in the triplet case. 
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To summarize, we have considered phases reahzed in microscopic weak coupling theory in 
multicomponent superconductors. We found that, for the most part, weak couphng mechanism 
in superconductors selects the most interesting magnetic phase, which is similar to A-phase 
in superfluid ^He. In this pairing state the time-reversal symmetry is broken. The properties 
of such superconductors are well known [7 10]. Similar to lm\k fcrromagnets, a pairing state 
which breaks time-reversal symmetry tends to form domain structure, with magnetic moment 
arising from circulating electric current in the domain walls. Thus, there is no net magnetic 
moment in the bulk. Perhaps, the most interesting phenomena for multicomponent super- 
conductors arise in magnetic field. An obvious consequence of magnetic field is the change of 
phase for a certain direction, which has to do with the symmetry [9,10]. This gives rise to very 
unusual topological excitations and vortex lattices [9,21,22], and an H — T phase diagram, 
which reflects the possibility of many different phases [6, 23] or an abrupt change of slope for 
Hc2{T) [16]. Another detectable property is the splitting of Tg under uniaxial stress [24]. 

Our results are summarized in table I, where we also list specific heat jump at Tc for simple 
Fermi surfaces, such as cylindrical or spherical, written as a ratio to s-wave BCS theory 
specific heat jump, (AC/C) = 12/7^(3) ~ 1.42. Our results for singlet superconductors 
are consistent with Ref. [19], where the BCS mechanism for multi-pocket Fermi surfaces 
was considered. The main conclusion of Ref. [19] was that it is the peculiar shape of the 
Fermi surface, such as its multi-pocketed structure, which determines the phase realized in 
BCS approach for a multicomponent order parameter and leads to a magnetic phase. We 
found that this statement is only partially correct. For all degenerate representations of the 
cubic group, except for the 2D representation Eg, the superconducting phase depends on 
the shape of the Fermi surface. The weak coupling approach, however, introduces certain 
constraints on the GL coefficients and possible phases. For the 2D representations of the 
hexagonal and tetragonal groups and the representation Eg of the cubic group, the magnetic 
superconducting phase in the weak coupling approach is determined by symmetry alone. 
The shape of the Fermi surface determines the specific heat jump at Tc and even the form 
of the gap function, since the basis functions for irreducible representations of point groups 
are not unique, but not the superconducting phase. The only requirement for magnetic phase 
to appear is that the superconducting order parameter transforms according to one of these 
irreducible representations. 

This work was supported by TAML at the University of Tennessee. 
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irr. 
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(rii,m) or 


Magnetic 
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group 


repr. 


class 


(rii,ri2,r}3) 


properties 


spher. 
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0(02) 
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DeiE) 
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FM 


0.7 BCS 


BCS 


Deh 


Elu 


De{E) 


(1,0 


FM 


5/6 BCS 


BCS 




E2g 




(1,0) 


FM 


0.7 BCS 


BCS 


Hexagonal 


E2u 


D2®R 


(1,1) 


non-magn. 


5/6 BCS 


BCS 






02(02) (SR 


(1,-1) 


non-magn. 


5/6 BCS 
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Eo 


Di(E) 
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0.7 BCS 


BCS 


Tetragonal 
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FM 


5/6 BCS 


BCS 
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